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Abstract 

The influence of the strong laser-driven vacuum on a propagating electromagnetic probe wave has 
been studied in detail. We investigate two scenarios comprising a focused probe laser beam passing 
through a region of vacuum polarised by an ultra-intense laser field. By splitting this strong field 
into two, separated, monochromatic Gaussian pulses counter-propagating in a plane perpendicular 
to the probe field axis, we demonstrate a leading order light-by-light diffraction effect that generates 
an interference pattern reminiscent of the classic double-slit experiment. We calculate the total 
number of probe photons diffracted as well as the number diffracted into regions where the vacuum 
polarisation signal is higher than the probe background. In addition, we calculate the induced 
ellipticity and polarisation rotation in the probe beam and show how, in the realistic situation in 
which the centres of the two strong fields are not exactly aligned, certain ranges of beam separation 
and observation distance may actually lead to an increase over the idealised case of a single strong 
laser beam. 
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I. INTRODUCTION 



That strong electromagnetic fields can modify the dielectric properties of the quantum 
vacuum has been known since the pioneering work of Heisenberg and Euler, Weisskopf and 
Sauter Quantum electrodynamics (QED) predicts that when the electromagnetic field 

strength nears the critical "Schwinger limit" required to spontaneously create an electron- 
positron pair with an electron of charge — e < and mass m within the reduced Compton 
wavelength A c = H/mc, a range of nonlinear vacuum polarisation effects (VPEs) should be- 
come observable. The corresponding electric field of -E'er = y/4jrm 2 c 3 /he = 1.3 x 10 16 Vcm -1 
in Lorentz-Heaviside units, would certainly allow one to access attenuative VPE processes, 
namely involving real electron-positron pair creation, whose rates would become large 
enough to be easily observed. Moreover, recent calculations show that these effects in the 
presence of loan fields can already be clearly observed at intensities orders of magnitude 
below critical values 4j-|7J]. Such VPE processes have also been discussed as a probe for new 
fundamental physics, with current limits clarified and new experiments proposed js jj]. 
Since these processes are exponentially suppressed, for the case in earth-bound laboratories 
where the electric fields involved are much less than E cr , it is refractive VPE processes 
involving virtual electron-positron pairs that are most likely to be observed. Photon-photon 



scattering is one example of a refractive VPE 13], which has already been carried out 
as Delbriick scattering involving virtual photons in the Coulomb field of a heavy nucleus 
where the atomic number Z is < 1/a and a = e 2 /AT(hc ~ 1/137 is the fine-structure 
constant, but has since eluded detection for purely real photons [15] . This effect could be 



measured by virtue of polarisation-dependent emission in four-wave mixing 



161 ] : by using 



the transverse-electric modes of plane waves to generate a resonant coupling in a waveguide 



171 ] or using sufficiently intense lasers to compensate for the small cross-section and 
inducing a phase-shift in lasers passing through one another fisl . jlo| (for a review of the ap- 



plications of relativistic lasers see 2CH22|). Other nonlinear vacuum effects include "photon 



acceleration" [23|; photon splitting in atomic fields 24j], pair plasmas 25] and laser fields 



261 ] ; as well as the corresponding reverse process of vacuum high-order harmonic generation 



e.g. in various laser set-ups 
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or in a mixed Coulomb and laser field set-up [29 ]. 
The corresponding critical bound of the magnetic field, B CT = V^rm 2 c 3 / he = 4.4 x 10 13 G 
can be surpassed by ultramagnetised neutron stars or "magnetars" 30[ (for a review on 
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X-ray pulsars, see 3JJ), which provide an inhomogeneous trigger for nonlinear effects such 



as vacuum birefringence and photon ray-bending [32j as well as photon-splitting 33|, |34 |. 
The possibility of laboratory-based experiments that measure second-harmonic generation 
within a constant inhomogeneous magnetic field, have also been considered [35|. The 
current PVLAS (Polarizzazione del Vuoto con Laser) experiment uses a slowly-varying 
magnetic field to attempt to detect refractive-regime vacuum-induced birefringence and 
dichroism through rotation in the polarisation of a probe laser wave [3^]. In addition to in 
a magnetic field, birefringence can also be induced in the vacuum by e.g. a laser field in 



this regime 
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40(. This latter scenario, and that of vacuum-induced diffraction are two 



examples of refractive VPEs which we further develop in the current paper. 



At the time of writing, the record for the highest intensity laser ever produced is held by 
the HERCULES laser and stands at 2 x 10 22 Wcm -2 4l|, seven orders of magnitude removed 
from the Schwinger limit intensity of I cr = cE^/2 = 2.3 x 10 29 Wcm -2 . We foresee that 
with the next generation of lasers currently being built, we will soon be in a much better 
position to test vacuum effects and so work with the quoted values for intended intensity 
and photon energy ranges in the coming decade. Examples of strong-field lasers are the ELI 
(Extreme Light Infrastructure) and HiPER (High Power laser Energy Research) facilities 
with target intensity values of 10 26 Wcm -2 42|, |43j. The PFS (Petawatt Field Synthesiser, 
44j), whilst planning a lower intensity of 10 22 Wcm -2 , will have a repetition rate of 10 Hz 
which could be more favourable in certain situations and provides an example of a cutting 
edge system to be commissioned later this year. Free-Electron Lasers (FELs) where 
undulating electrons provide the lase medium, such as the XFEL (X-Ray Free-Electron 
Laser) and the LCLS (Linac Coherent Light Source) could also be used to polarise the 
quantum vacuum, especially with the so-called "goal" parameters quoted in |45|. However, 
a further application of the FELs, one which could be reached sooner, would be as probe 
field lasers, whose alteration when passing through vacuum polarised regions could be 
measured. The XFEL and LCLS would be ideal for measuring refractive effects, which are 
in general proportional to laser frequency, as they allow continuous adjustment of the probe 



wavelength down to a minimum of 0.1 nm and 0.15 nm respectively 46|, |47|. 



This paper concerns itself with laser-induced vacuum-polarisation effects in the spirit of 
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40| . Here, the change in polarisation and ellipticity of a planar Gaussian probe field passing 
through a region of the vacuum polarised by a perpendicular standing wave formed by two 
counter-propagating, ultra-intense (J > 10 23 Wcm -2 ) Gaussian beams, was calculated up 
until the point where probe defocusing becomes important. We compare and expand upon 
this simple set-up with the following enhancements: 

i. The two counter-propagating strong field wave triggers for VPEs are separated in the 
plane perpendicular to their propagation, modelling a more realistic situation. This 
makes sense first from an experimental point of view, to know how VPEs are sensitive 
to laser alignment, and second allows us to derive an interference effect as different 
parts of the probe beam pass through the "double-slit" -like, vacuum-polarised region. 
This will even turn out to increase probe-beam polarisation rotation and ellipticity. 

ii. Defocusing terms were introduced into the probe beam and a corresponding update 
to the vacuum-induced ellipticity and rotation of the probe polarisation. This extends 



the limited range of detector distances in 40( where we could have compared theory 



to experiment, as our new expressions are also valid in the far-field diffraction zone, 
where they converge to a non-zero value. 

iii. The electric field generated in the new set-up by the current of the polarised vacuum, 
which we henceforth label the "diffracted field," was also calculated in the probe 
beam's transverse plane, allowing us to again model the more realistic situation where 
a detector is placed off-axis, in regions where the diffracted field, having a wider spread, 
is larger than the probe background. 

Throughout, we will make the analogy with the single- and double- slit diffraction 
experiment. According to Babinet's principle, the diffraction pattern generated by light 
passing an opaque obstacle is the same as that for light traversing an aperture with the 
same shape as the obstacle j^sj]. Regions of the vacuum polarised by the two strong-field 
laser beams then represent "translucent" obstacles for photons in the probe beam, having 
as we will show, a non-trivial polarisation and magnetisation. Unlike the typically sharp 
two-dimensional slits used in demonstration experiments, the strong lasers, being Gaussian 
in beam profile, form smooth, three-dimensional slits. One consequence of this will be that 
no single-slit fringes occur in the far field. However, the probe photons scattered from each 
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strong beam will interfere with one another, and in this way, we will have a "double-slit". 
As the scattering of probe photons occurs with such a small probability, the complete 
double-slit pattern will only be observable when the background of probe photons passing 
unperturbed through the apparatus, is subtracted. At the detector, the total field will in 
general consist of the probe signal plus the vacuum contribution of the scattered photons. In 
calculating the interference between these two fields, we demonstrate a new diffractive effect 
of the polarised vacuum, which we accentuate by forming the double-slit-like experimental 
set-up. In addition, we also compare polarisation results with a second beam geometry, 
namely that of the probe propagating anti-parallel to the strong field, which we label the 
"double-shaft" (or "single-shaft") set-up. Our results are complementary to findings in (37) ]. 
which focus on only the pure diffracted intensity for a different laser geometry. 

The paper is organised as follows: in section II we first introduce the Euler-Heisenberg theory 
upon which the results are based and the range of experimental scenarios we consider; then 
follows in section III an analysis of the first part of the results, the intensity of the bare 
diffracted field and the time-averaged difference in 'probe + diffracted' signal in terms of 
number of photons; the second part of the results deals with the change in rotation and 
ellipticity of polarisation for both of the two geometries and the paper is concluded with a 
recapitulation of the main results. 

II. THEORETICAL BASIS 

A. Leading-order vacuum current 

By making two basic assumptions, we can drastically simplify the interaction terms occurring 
in our field theory [49j. From the assumption that the photon energies involved are much 
lower than the electron rest energy, follows that the loop contribution of spatio-temporal 
extent can be consistently regarded as a single local interaction point, thereby allowing us 
to write down the so-called Euler-Heisenberg local, point Lagrangian density, C, which we 
will use to describe vacuum polarisation effects. Secondly, from the aforementioned laser 
intensities which are either currently available or scheduled for the future, we will work 
with the comfortable assumption that field strengths are much lower than critical values. 
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This then allows us to use the weak-field expansion of the Euler-Heisenberg Lagrangian, the 
leading order of which (in a system of units adopted henceforth, h — c—1) reads: 

£ = \(E 2 - B 2 ) + [(E 2 - B 2 ) 2 + 7(E • B) 2 ] , (II.l) 

for electric and magnetic fields E and B and their square moduli E 2 = E • E, B 2 = B • B 
respectively. Extremising the action with respect to the vector potential corresponding to 
these fields, we achieve the following wave equations of motion, for an induced vacuum 
current J vac : 

V 2 E-<9 2 E = J vac = V A d t ~M — V(V • P) + <9 2 P, (II.2) 

4a 2 

p : = ^— \2(E 2 - B 2 )E + 7(E • B)B , (II.3) 

45m 4 

M^-^1[2(^-^)B-7(E.B)E]. (II.4) 

There are many similarities one can draw between birefringent solid-state materials and 
the behaviour of the vacuum under intense electromagnetic fields. Direct from the above 
wave equation Eq. (111.21) . we can liken the vacuum current to one representing the response 
from such a birefringent material, that is to say, labelling P its polarisation, and M its 



magnetisation 



551 ] . As these are functions of both E and B, the inhomogeneity in our fields 



which here plays a central role, is included at this point. 



Definition of experimental scenario 

In this paper, we will focus mainly on the double-slit set-up sketched in Fig. (TJ and include 
only a summary of the polarisation results for the "single-shaft" set-up corresponding to 
a head-on probe and strong-field collision, towards the end of the article (we use here the 

n 

label "shaft" contrary to in [37J, where it was also labelled a "slit"). For our double-slit 
set-up, two tightly-focused (we assume the diffraction limit has been reached i.e. focused 
down to a wavelength), counter-propagating monochromatic strong field Gaussian pulses 
polarising the vacuum, with waists w ,o centred at (x, y) = (a, b) and (x, y) = (—a, —6), elec- 
tric fields E ,i(x, y, z, t), E 0i 2(^, y, z, t) and wavevectors k = (0, 0, — uq) and — k = (0, 0, uq) 
respectively, are permeated simultaneously by a broader and weaker linearly-polarised, trans- 
verse monochromatic Gaussian probe field, polarised at an angle 9 to the x-axis with waist 
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w Pi o, electric field E p and wavevector k p = (0, tu p ,0), which will gain a diffracted compo- 
nent, whose intensity and polarisation will be measured some distance, r& away. Using the 
Gaussian beam solution from 201 ]. in the effective interaction region we therefore have the 



following: 

E (x,y,z,t):= E 0jl (x,y, z,t) + E 0i2 (x,y, z,t) x, (II.5) 

UqZ (x — a) 2 + (y — b) 



Z* + Z- 



E 0il (x,y,z,t) :=E 0fi (x-a,y-b,z) sirM + u t + u z - <j> &0 (z) + 

_, / . _, . is./, ,/\ w 2 (x + a) 2 + (y + 6) 2 \ 
#0,2 (a, y, 2, *) := ^o,o (z + a, y + b,z) sml^o + - w 2; + g ,o(z) „ ^— i J 

P + w p t - w p y + g , p (y) J — 2 2 J (x cos 6» + z sin 0) 



(II.6) 

where we have defined respectively the strong and probe fields amplitudes Eo,o(x,y, z), 

E p fi(x, y, z) with their maximum values E /y/2 and E p <C E , as: 

E Q e -(x 2 +y 2 )hl e -(^ 2 +z 2 )/wl 
E 0fi {x, y, z) := — ————— t E pfi (x, y, z) := E p 



v/2v/l + (zM,o) 2 ' P ' UV Vi + Cf/fr.p) 2 ' 

The square of the waist of focusing is defined from beam parameters as Wq : = 
w o. 

(l + (z/z rt0 ) 2 ), w 2 :— w 2 (l + (y/i/r,p) 2 ), where w P;0 > w ,o, with Rayleigh lengths defined 
in the usual way, z T ^ Q = uqWq Q /2 = 7rw 0i o (as we have assumed w 0; o = Ao), y T , P = w p w p ,o/2, 
and the Gouy phases respectively <j>g,o(z) = taxr 1 (z/z Ty o), 4> g , P {y) = tan _1 (|//?/ r]P ). The 
fields in Eq. (1II.5j) are chosen as a first-order approximation to the solution of Maxwell's 
equations in vacuum (see e.g. [20J, p.p. 64-65 for details on the higher order terms in this 
expansion), being an expansion in the small parameters e z = w 0l o/^r,o — Ao/(vrw 00 ) rj 
and e y = w Pi0 /y r ,p = A p /(7rw P)0 ) <C 1 (as by definition the probe is not intensely focused 
and so w p A p ). Therefore, throughout this calculation, we are working to an accuracy 
given by the largest term neglected in the expansion, e z . 

The magnetic fields consistent with this level of approximation are then: 

B (x,y,z,t) = -[E 0tl (x,y,z,t) - E 0)2 (x,y, z, t)]y, (11.7) 

B p (x,y,z,t) = E Pj o(x, y, z)(xsin 9 — z cos 6) x 

/ oo y x 2 -\- z 2 \ 

sin ip p + u p t - u p y + <f) gtP {y) — =- — =- . (II.8) 

V 2 y 2 + yl p J 



Since the probe field's strength is much lower than the strong field's, we regard terms 
oc Ep,Bp and smaller as being negligible in Eqs. flll.3j) and flll.4j) . In addition, as we are 



only interested in the effects on the probe field, we drop terms which only depend on the 
strong field oc E^,Bq. With E = E + E p and B = B + B p , the vacuum polarisation and 
magnetisation then becomes: 
4a 2 



P 

M 



45m 4 
4a 2 

45m 4 



2(E ■ E - B • B )E p + 4(E • E P )E + 7(E ■ B p )B j , (II.9) 
2(E • E - B • B )B P + 4(E ■ E P )B - 7(E ■ B P )E 1 . (11.10) 



Diffracted field off-axis 

We will focus on the diffracted electric field E d (r d ,t), generated by the polarised vacuum 
current in Eq. (1II.2j) . at a displacement r d , from the centre of the interaction volume (centred 
at the origin of the co-ordinate system), in the direction of propagation of probe beam. Using 
Green's functions to solve the inhomogeneous wave equation driven by a current J(r, t), we 
have, in general, E d (r d) i) = — 1/(4-71") J d 3 r|r d — r| _1 J(r, t — |r d — r|) j^]. It can be seen 
from the definition of our current in Eq. ( 111.21) . that we are going to have terms cubic in the 
electromagnetic field, which means cross-terms between our probe and strong fields in the 
interaction region. As our waves are monochromatic, we see that photons of discrete energies 
Up and Up ±2o;o, etc. will be produced. Photons with the latter energies are evanescent and 
can therefore be neglected, which turns out to be equivalent to averaging the expression in 
time. We Fourier transform our current in time in order to use this discreteness and then, 
as we are only interested in effects in the probe, will later set the frequency u to u p : 



E d (r d ,u) = ~ f d 3 r dt\v A d t M - V(V • P) + 9 2 p1 [ , r|+t)J . (11.11) 

4vr J L J |r d — r| 



It will be useful to expand the exponential using the assumption that the detector is 

placed much further away than the dimensions of the interaction volume, taken as the 

standard deviation width of the beams. Using wo,o < w Pi o <C r d , and then assuming 

(wo,oM P )(w Pj o/rd) 2 , (w Pi o/A p )(w Pi o/Vd) 3 <C 1, we can curtail the expansion to: 

exp I" — iu(\r d — r| + 1)] 1 \ / r 1 , l91 \ | 

-JL1 — vi g_ i — n w — eX p J _ i0J t rd _ fd . r+ | fdAr |2 +t \ I (IL12 ) 

|r d — r| r d I V 2r d / | 

By retaining the quadratic co-ordinate terms, we indicate that we'll be working in the Fresnel 



regime. We can then split Eq. (III. lip into three integrals and integrate by parts to remove 
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surface terms, which we assume, using Gaussian expressions, tend to zero at the boundaries. 
This leaves us with: 

,2 . 



E d (r d , u) 



' exp [ — iur d \ 
Aura 



d 3 r dt (M A f d + P — P ■ f d f d ) x 

1 



(11.13) 



exp 



iw(f d ■ r - — |f d A r| 



We expect the main vectorial contribution to the probe from the vacuum polarisation and 
magnetisation to be in the x and z directions, i.e. the directions of the probe and strong 
electromagnetic fields. When we substitute our particular scenario using Eq. (III. 91) and Eq. 
(III. 101) into the above equation and then Fourier transform back into (r d , t) co-ordinates, we 
achieve the following: 



P f r A F'fr ^xp[i(-u p r d + u p t + ip p )] exp[-i(-u p r d + uj p t + ip p )) 



E d (r d ) := 



h aE p 
I cr 45A2r d 



2i v ' 2i 

(V! + V 2 )ui + (V 3 - V 4 )u 2 + ( H J u 3 



(11.14) 



i=l 



where the volumes Vk, and the vectors, Uj are defined as the following: 



d 3 r exp 



+ iui 



x 2 + y 2 + z 2 xx d + yy d + zz d 



2r d 

(xx d + yy d + zz d f 
2r 3 



+ y) 



r d 
x 2 + z 2 



(11.15) 



w 



p.O 



Zl 


:= exp 


- 2 

- w o 




:= exp 


■ *0 




:= exp 


- 2 

■ ™l 




:= exp 


- 2 
■ w o 



exp 
exp 



- 2i\u Q z - <j) g fl(z) + 
2i(cu z - 4> g ,o(z) + 



i+(z/ Zr r 

u z(x 2 +y 2 + a 2 + b 2 ] 
2{z 2 + z 2 ) 
oo z(x 2 + y 2 + a 2 + b 2 ) 
2(z 2 + z 2 ) 



(x - a) 2 + (y- bf 



-zUx + aY + iy + b)' 
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/ 



Ui 



(1 



'Ha' 
'd - 



cost 



^(^-cosO + ^sin( 

r d V r d r d 



— sin 6 -\- — cos 9 — cos# + M sin( 

r d r d r d V r d r d 



''a - 



sinO - Sd(^L cos6 + am 



u 2 



V 

^ ^cos0 + ^^sin0 ^ 



''d 



r d ' 4 r d r d 

It 

4 V \ r d • 

-^cosfl + l^sinfl 

r d 4 r d r d 



U 3 



(l-(f^) 2 )cosfl 

4 r d r d r d 

v |^sin0-^cos0 

\ 4 r d r d r d 



/ 



The main contribution from the integrals, Vj, will be within the widths of our laser beams 
and so we can regard x, y < w , z < w P)0 . If we evaluate the expression at these values, the 



probe amplitude defocusing terms become (1 + {y/y T ^ 



1 - (l/2)(y/?/ riP ) 2 , and we see 



that the correction (l/2)(y/y rp ) 2 <C e z (the accuracy of our computation), and is therefore 
negligible. Moreover, considering the defocusing phase terms, when we assume these ranges 
for x,y and z throughout the integration, <j> g)P (y) ~ (y/y r ) 1 an d the final defocusing 
term u p y(x 2 + z 2 )/[2(y 2 + y 2 p )} ~ 2A p wo,o/(ti"w 2 )0 ) <C 1 for the realistic parameters that we 
take for our lasers, defined later on Sec. IHIl P. [121 Therefore, to be consistent with our 
beam expansion, we have considered all probe defocusing terms as constant (w p (y) — >■ w Pj0 ) 
within the integral above Eq. flII.14j) . Whenever the probe occurs explicitly in expressions 
outside the integrals, the full space-dependence will be used. 



From this integral Eq. ( III. 151) . it can be seen that the x and the y co-ordinates can 
be integrated out to give just an integral in z (see Eq. 11A.1I) in the appendix). On 
inspection, we notice certain factors in the complex exponential of the integrand constrain 
the diffracted field to be sharply peaked around Xd/^d ~ and yd/r<i ~ 1, agreeing with 
physical intuition. Taking the limits — > yd = (0,y^,0); a, b — > 0, we can easily recover 
the expression for an on-axis measurement of a single strong beam + probe collision, 
given in [40] . Similarly, we can derive the diffracted magnetic field Bd(ra,u; p ), and using 
Maxwell's inhomogeneous equations again, show Bd(rd,Co> p ) = k p A Ed(rd,Co> p ) to within 
our calculational accuracy and i.e. that the flow of energy described by Poynting's vector 
goes as Sd(r d ,w p ) = E d (r d ,(X> p ) A B d (r d , w p )/2 = |E d (r d , w p )| 2 k p /2, which simplifies our 
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calculation of the intensity pattern. We have assumed our earlier conditions on w ,o and 
w Pj0 that we used in Eq. flII.12j) . as well as that (2/d/fd) 2 ~ 1- 



One can question how sensitive these results are to being able to align the strong-field lasers 
parallel to one-another, by considering them being focused from a distance away by two 
large mirrors. For a small rotation 6(f) of ko,! and ko,2 around the x-axis in the directions 
±y respectively, one can show for W 01 — B 0i iCos5</> + E 0jl smScf)z, ^02 = B Oi2 cos50 — 
£0,2 sin 50 z (E' ,i = E o,i,Eo i2 = E , 2 ): 

PW) = P-¥^(^o 2 i-^o 2 2 )cos^ y + O[(50) 2 ], (11.16) 
lfio- 2 

M(54>) = M + ^^^-^^cos^y + OK^) 2 ]. (11.17) 

Keeping within the aforementioned bounds in the detector-plane, the corrections in 5<p 
cancel in the combination (M A fa + P — P?d ?d) meaning corrections to Ed ~ (<50) 2 , 
which implies, envisaging 5<p ~ 0.1, that the parallel idealisation is sufficient to within the 
accuracy of the present treatment, 

A further consideration would be what role the finite length of the strong-field beams plays 
during the passage of the probe beam. Both the diffracted intensity and polarisation effects 
that we will study are proportional to the intensity of the strong-field and therefore the 
corresponding longitudinal distribution is ~ 1/(1 + (z/z Ti o) 2 ). In the absence of a well- 
defined decay length, we take the effective length to be that at which the intensity falls to 
below a tenth of its initial value, giving an effective length of Iq = 3z ri o = 7.5 /mi. A finite 
pulse length leads to consideration of the temporal envelope. For the case of a Gaussian 
beam, the leading temporal correction should be of the order 1 / ooqTo which we have already 
specified, through our assumption of monochromaticity, to be <C 1. If we ensure that the 
strong-field pulse-length r is such that cr > 2l then the deviation should be negligible to 
within our level of accuracy. We therefore choose ctq w 2 x 21q, with tq = 100 fs, which will 
limit the maximum strong-field intensity obeying io r o^4 = £ for a fixed laser energy £ and 
focus area A. 
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III. ANALYSIS OF THE RESULTS 



We present results that follow from the numerical evaluation in matlab of the one- 
dimensional integral for Ed(rd) given in Eq. (jA.ip in the appendix. 



The results are presented in two sub-sections for i) intensity and ii) polarisation. These are 
further divided into the form of intensity along the x-axis, along the z-axis and in the x-z 
plane; following which we explain the polarisation rotation and ellipticity expressions along 
the probe propagation axis. 



Intensity measurements off-axis 

The nonlinearity of the vacuum brought about by the two strong-field waves generates the 
diffraction patterns one would expect from a refractive solid-state material. The integral 
expression of Ed(rd, oS) in Eq. fill. 131) allows one to interpret the effect at hand as an example 
of Fresnel diffraction, including as it does, squared co-ordinate terms in the exponential. 
Satisfying the inequality: (w ,o-2r,o/^p r d)(a ; d + ;2; d)/' r d C 1, we can neglect the xz cross terms 
in the exponential and form two independent diffraction parameters, £ x = WQ /A p rd and 
£ 2 = Wp /A p r d , for which £ X)Z ^> 1 implies taking the near-field limit, while £ X)Z <C 1 implies 
we can take the far-field limit and hence the Fourier transform of the transmission function 
48]. 



Analysis of field diffracted onto the x-axis 

Numerical evaluation of the leading-order QED contributions to the field diffracted onto the 
x-axis is shown in Fig. [2j For our probe beam, we take the X-FEL at DESY in Hamburg, 
for which we have 80 GW in a pulse of length 100 fs, of 0.4 nm wavelength radiation focused 
into a waist w Pi o = 100 /xm 4^]. We maximise the intensity of the diffracted field by 
setting 6 = n/2. In addition, we take for our strong- field beams parameters from the 
upcoming ELI and HiPER facilities i.e. Ao = 0.8 /zm and assume that they can be focused 
up to the diffraction limit i.e. that Ao = wo,o (although the consequences of focusing to a 
width of a few wavelengths are not drastic for the results). As already discussed, we then 
choose a pulse duration of r = 100 fs to satisfy cr ~ 4Z and a total peak intensity of 
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/ = 10 24 Wcm- 2 . The strong beams are separated by a/wo,o = 6 and observations made at 
yd = 1 m. In Fig. |2j we plot the diffracted field intensity, which clearly shows a familiar 
squared cosine, with a symmetric, decaying envelope function, similar to the square of the 
Fourier transform of a double-slit transmission aperture. This result is expected if one 
notes that with the above numerical parameters, the diffraction parameter £ x along the x 
direction is much smaller than unity. We also note at this point, that separation of the 
strong beams in the direction transverse to detector plate has in general no observable 
effect on our numerical results for intensity, which can be understood intuitively. As the 
vacuum signal Ej is created in phase with the probe, the total phase difference between 
sources of vacuum waves separated by 2b in the longitudinal direction is 2bu p (l — y&/r<i), 
as can be seen from Eq. (111.120 in the far-field. Setting z& = for simplicity, the condition 
to be fulfilled for a corresponding first minimum would be A p /2 = b(xd/rd) 2 - Since we are 
using an X-ray probe, for realistic separation of the strong-field beams of the order of a few 
multiples of w ,o, the first minimum would occur for values of Xd far outside our detector 
plate. Moreover, for a finite strong-beam x-separation 2a, any additional separation of the 
beams in the y-direction, will neither create an appreciable separation perpendicular to the 
diffracted wave vector. For these reasons, we can disregard b and set it equal to zero in this 
section. 



In the present case x&, Zd <C yd and this implies that the terms proportional to the vectors Ui 
and u 2 in Ed(rd) are negligible. Also, we notice that for the typical situation, w ,o+a <C w P; o, 
the cosine term formed from the integrals Vi + V2 can be neglected when both: 



27TWp ;0 4 



2 - 2 



N \ 1 + Mr and « L (IIL1 ) 

A \ y Vd\ J 2A P r d 

These observations then considerably reduce our diffraction integral in Eq. (III.14P to just: 

Ed ( rd ) ~ TF^-rW+^W (ni.2) 

The full Fresnel-like diffraction integral which couples the x, y and z co-ordinates together, is 
unwieldy when attempting to garner qualitative information. Assuming £ x <g 1, the Fresnel 
integral will produce a diffraction pattern with the same shape as if we took the Fourier 
limit. In this way, by performing the integral in x in V3 and V4 we obtain for the diffracted 
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field intensity, Id(rd, u> p ) = |Sd(rd,w p )| that: 



^d(^d, J/d, = 0, u p ) ~ J Pi0 exp 



{x d /r d ) 



2-, 



COS 



u; p a (x d A"d) 



0". 



AjV2 

7TW n ' 



(III.3) 



with J p o = Ep/2, which is what one would expect from the Fourier transform of a Gaussian 
convoluted with two delta functions. The cosine term originates from the interference be- 
tween the vacuum current generated in the two slits, and the Gaussian is the effect of the 
single-slit shapes of both the strong beams. We use the ~ sign to emphasise the illustrative 
nature of our arguments, as although the fringe positions are correctly predicted, the single- 
slit shape is incorrect, as seen in Fig. El This is an example of a consequence of non-trivial 
beam geometry, for which the full three-dimensional integration must be performed. 



Analysis of field diffracted onto the z-axis 

An example of a diffraction pattern in the z-direction is shown in Fig. [3j The numerical 
parameters are those used in the above case but with a/w 0> o = and now the reverse 
situation x d = and z& <C y d . From this figure, we see that the intensity pattern is formed 
by a central peak, of width w 50 /mi, and two smaller exponential-shaped peaks some 
distance away. Concerning the central peak, when we consider that the amplitude of the 
strong field along the z-axis, and hence the "vacuum transmission aperture" is governed 
by the factor 1/(1 + (z/z r)0 ) 2 ), we see very clearly that the diffracted electromagnetic field 
does not result from the aperture's Fourier transform, which would have been a decaying 
exponential, symmetric about the origin, i.e. the wrong shape and with a smaller width of 
about 10 /im. The presence of the two peaks can be described by the diffraction-grating-like 
sinusoid along the z-axis. That the simple Fourier analysis applied in the previous case does 
not work here, is already clear from the diffraction parameter £ z ~ 25 not being smaller 
than unity. 



We wish to again explain our diffracted field qualitatively, but now how it, and so how 
Eq. ( III. 141) . depends upon the z co-ordinate in the detector plane, z&. We can see from 
Eq. (jA.ip how the decay of the integrand in the z-direction is controlled by the softcore 
1/(1 + (z/z rj o) 2 ) term. The importance and presence of this term prevents us separating out 
the z from the x and y integration variables and hence the Zd from the x d and y d detector- 
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plane co-ordinates in the integrands Vk, as we managed to do in the previous case. However, 
if we consider that Xd = and that Zd <C yd, we can again neglect in Eq. (III.13j) the terms 
proportional to the vectors ui and u 2 . Unlike the previous situation however, the condition 
Eq. ( IHI.ip to neglect the integrals Vi and V2 is not fulfilled for arbitrary z& <C ?"d and they 
are accordingly not negligible. It can be shown by performing an analysis similar to the one 
in the previous case that the integrals V3 and V4 give rise to the central peak (with width 
w Pi o/2 = 50 fim) while the integrals Vi and V2 give rise to the secondary smaller peaks 
located at Zd = =F2rdA p /Ao = =Fl000/im. Therefore, the secondary peaks originate from 
the standing wave of the strong field, which the probe experiences as if it were a diffraction 
grating, and is another example of the effect of non-trivial beam shape. Similar arguments 
leading to Eq. ( jIH.3p . retaining the quadratic terms in the exponential give the dashed line 
in Fig. [3] and again show good agreement. 



Single-slit pattern 

We have now seen from some results that a consequence of the non-trivial strong-field beam 
shape is a deviation from the ideal double-slit analogy. As mentioned in the introduction, 
this mainly affects the interpretation of each strong- field laser as a single slit. We can 
illustrate the difference brought about by smooth edges when we consider diffraction from 
a single slit of dimension 2l x by 21 z , centred at the origin. The diffracted electric field in the 
far-field, -Ed,Rect can be calculated via Fourier transformation of the aperture function: 

Ed Rect oc / dx dzexp [ —iu—x — iu—z ) Rect ( — ) Rect [ — ) , (IH-4) 
'-00 \ r d r d J \l x J \IJ 



oc 



where Rect(x/a) equals unity only in the region x g] — a,a[, being otherwise zero. This 
gives a diffracted intensity I d ,Rect oc l-^d ,Rect I • 



sin 2 (o;a: d /r d ) sm 2 (ujz d /r d ) 
(uxd/rd) 2 (uz d /r d y 



am i^u/Ad/ ' d) 0111 \yj^dl 1 dj /TTT r i 

/d,Rect OC j_ - j- ^ — j- , (III.5) 



which gives the familiar single-slit minimum conditions (n+l/2)X = 2l x Xd. n /^di (h+1/2)A = 
2^^d,n/ r d, for n G Z and A = 2-k/uo. For our "Gaussian" slits, our diffracted electric field 
^d.Gauss becomes: 

£d,Gauss oc / dx I dzexp (-iu—x - iu—z\ exp ( 2 * 2 ) j 2 - (III. 6) 

J-00 J-00 \ r d r d / V Waov- 1 - + z )/ 1 + z 
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This can be analytically evaluated after setting z d = 0, giving an intensity: 

-fd,Gauss oc exp 



V 2r d J 



K 2 



1 / o;wo,oa: d \ ; 

2 V 2r d / 



(III.7) 



where i^o is the zeroth-order modified Bessel function of the first kind and is monotonically 
decreasing, i.e. without fringe structure. As other terms introduce only a finer structure 
and as the final integration in y would also be over a smooth function, we see that no pe- 
riodicity arises from our single-slit diffraction pattern, which is consistent with numerical 
results. Beyond the far-field limit however, an interference-like deviation would be expected 



to develop. One example of this was calculated in [5l|, where a relativistic Gaussian elec- 
tron wave-packet in the Coulomb field of some highly-charged ions acquires an interference 
pattern structure when placed in an intense laser field. 



Resultant intensity difference off-axis 

For the relevance to experiment however, instead of just plain diffraction theory, we will 
be more interested in studying the difference brought about by vacuum polarisation effects. 
With () denoting an average over a laser cycle, the difference can be shown to be: 



^tot — Ip 



Ipd 

h 



((|E P + E d | 2 ) 

Ipd + -^d, 

Jo «/ p , ex P 



(|E P | 2 )) 



-(^d + ^d)/ w 



I a 180ttA2 
(|£ d | 2 >; 



2/d\A + (yd/y r y 



-(V l sin?7— V r cos77)-(xcos0+zsin0) , (III. 8) 

(III.9) 



V r := real(V), V := imag(V), 
V = (Vi + V 2 ) Ul + (V 3 - V 4 )u 2 + ( J2 V< ) u 3 , 



ri = tan 1 ( — ) — 

\VvJ 



i=l 



(111.10) 

(III.ll) 

(111.12) 



2 y 2 T +y 2 A 

We can evaluate this expression on an x d -2 d grid at a fixed distance y d , and calculate the 
differences in photon rates brought about by the polarised vacuum. Our procedure was 
to make many such grids that became ever-finer, so that we could see how the integral and 
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i.e. the predicted number of photons per shot converged. Since the diffraction pattern must 
be smooth, this number should then yield a reliable value. We took the same following 
parameters of a typical experimental run: y d = lm, a/w = 12,6/w 0i o = 0,w 00 = A = 
0.8 /mi, Iq = 10 24 Wcm -2 , giving the patterns shown in Fig. (the parameters of the probe 
field were those already employed in the previous examples). We focus on the diffracted 
photons described by the interference term between E p (rd, t) and Ed(rd,£). This term 
spreads out in the x-z plane with a width a/2 larger than that for the probe field, as the 
multiplying diffraction signal has a much wider overall decay, so there exist regions in which 
the ratio of diffracted to probe signal is favourable which can be seen on the log-plot of 
the total difference due to vacuum signal over the probe background Fig. lib. At the same 
time, moving too far from the centre of the pattern will reduce the intensity to the point 
where nothing can be detected. If we consider drilling a hole of radius p into the centre of 
the detector and approximate the decay of J p d to come entirely from the probe Gaussian, 
considering the single strong-beam scenario in order to maximise signal, we can obtain limits 
on p: 



iV, 



pel 



1/2 < < 



lniV, 



1/2 

(HI.13) 



for total incident probe and cross-term diffracted photons N p , N pd . This agrees with the 
intuitive notion that to stand any chance of measurement, the signal must be larger than 
statistical noise^ from the background, which if modelled with Poisson statistics implies 



561 ] . We can either fulfil this condition that the vacuum signal is larger than 
the minimum background noise over the entire plate, or we can consider measuring counts 
only in regions where N p d(p) > N p (p). In both cases, the number of diffracted photons 
will simply increase with probe intensity, whereas as N p depends only on the probe laser 
energy and wavelength A p and so for larger probe intensity, we can easier fulfil both bounds 



on p in Eq. ( 1III.13|) . First setting y$ = 50cm and the still at ELI comfortably attainable 



Iq = 5 x 10 24 Wcm -2 , for a probe focal width of 8 pm, we achieve iV p d = 7.5 x 10 7 from 
N p = 8.0 x 10 12 probe photons per shot. Secondly, we can plot how N pd (p) varies with 
hole radius, and for a tighter probe beam focal width of w Pi o = 3.6 pm, which, in the 



light of recent results of focusing hard 20 keV photons to a width of 7 nm [52|, we expect 
to be attainable in the near future, we achieve the dependency shown in Fig. In the 
region iV p d(p) > N p (p), taking into account the efficiency of commercially- available CCDs 

17 



for A p = 0.4 nm or 3.1 keV photons (> 90% [531]). we expect approximately two diffracted 
photons to be measurable per shot of the probe beam. This can then be compared with 
results for J d which, by not being subject to the probe Gaussian envelope, has a much wider 
spread, and is possibly easier to measure as reported in j^], with the caveat that an optical 
probe beam was used with a total energy 2.5 x 10 3 larger than in the present X-ray case. 



Polarisation results (double-slit) 

This section concerns itself with the induced ellipticity and rotation of the probe polarisation 
due to VPEs, that can be measured on the probe beam's propagation axis. Setting Xd = 
Zd = 0, it can be shown that new expressions that incorporate defocusing terms in the probe, 
for the polarisation, ip and ellipticity, e, are given by: 

120A2 Jcr^W.P Vj ^ ' 

a sin 29 I A / V£ V, 



where in the limit of a, b, — > 0; y TtP — > oo, we again recover the expression in the original 



paper 



401 ] . We also note that the introduction of experimentally relevant defocusing terms 



in the probe, produces the more realistic and expected result that lim^^oo {■?/>, e} ^ {0, 0}. 



Varying yd with Xd/rd, Zd/rd = 0, and setting 9 = 7r/4 to maximise the effect of the polarised 
vacuum, we show a demonstrative plot in Fig. EK, for how ip and e vary for a fixed strong- 
field beam separation, a/wo,o — 12, with the other parameters the same as in the previous 
examples. The first difference we note is that in comparison with results from [40], for 
Vd "C y r , P , polarisation and ellipticity oscillate rapidly and there are sizeable ranges where 
both are larger than that for previously derived results. For the choice of parameters in the 
plot, y rj p ~ 80 m, and so if we keep within this range, i.e. disregard the effect of defocusing 
terms, we can clearly ascertain the improvement brought by separating the strong-field 
beams. This perhaps counterintuitive result can be shown to be consistent with our analysis 
by following through these conditions on the detector-plane co-ordinates and studying the 
form of the integrals 1/^4 which appear in our expressions for ip and e (we can once more 
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disregard the contribution of Vip)' 
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p,0 J 



-2(x =F a)' 



(111.16) 
(111.17) 

(111.18) 



From Eq. 1 T|) we can see more clearly, that under these conditions (especially as y r p 3> 
2(w 0j o + o)), since there is no other structure in the y-direction, b becomes an inconsequential 
parameter when measuring polarisation and ellipticity, just as it was for the diffracted field, 
and will likewise be set to zero. By separating strong-field beams in the x-direction, we see 
that we only produce an effect on the x- integrals, J Xi ±. When considering the contribution 
from the first complex exponential factor in Eq. (lIII.18p . for a fixed a, in varying y^, we 
vary the overlap this factor's real cosine and imaginary sine functions with the other two 
Gaussian integrand factors, which have maxima at x = and x = =Fa respectively. Hence 
some values of form local maxima in if> and e, and due to the trigonometric nature of 
the varying function, we have the oscillating shape in Fig. [6j However, in the limit y$ — > 0, 
(taking into account all V's), both of these values tend to constants: 

ay/n I w ,o f~ 2a2 \ ■ no mum 

We also show how ip and e depend upon beam-separation a in Fig. [TJ and can show 
consistency by using the same arguments as above for the dependence on y&. In varying a, 
the first two factors in Eq. ( 1III.18P act as fixed peaks, whereas the final Gaussian term is 
moved to place its peak x = =Fa, at such a position which could be used to maximise the 
integral. We recall from Eq. (IIII.14P and Eq. (IHI.15P that ip and e contain mixtures of both 
the real and imaginary part of this integral. When considering the contribution from the 
imaginary part of the integrand, we see that the maximum of the first complex exponential 
factor, i.e. of the sinusoidal, will not occur at the origin, unlike that of the second Gaussian 
factor, and hence in order to maximise this integral comprising three functions we should 
place the peak of the third function somewhere between the peaks of the first two, which 
corresponds to a value a ^ 0. Moreover, as the first sinusoidal factor is periodic, and has 
a wavelength much smaller than the width w Pi o, of the Gaussian which multiplies it, we 
should have a series of maxima in both ip(a) and e(a) which decay slowly with a (see Fig. 
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[7]). For the case y d <C y T , P , our explanation would on the one hand predict that the value of 
e(a) would initially rise as a increases, and on the other hand justify the maximum of il>(a) 
being very close to the origin, and hence that ip(a) would decrease as a initially increases 
from 0. These results can be further confirmed via differentiation under the integral in Eq. 
(IIII.16p . and are exactly what we observe in the numerical evaluation depicted in Fig. [7J 
This increase is another reflection of the role of Fresnel terms in a non-trivial beam geometry. 
From numerical analysis, the polarisation and ellipticity were found to increase by a factor 
of 1.4 over a = values. 



Polarisation results (single-shaft) 



We want to consider here briefly a different field configuration, in which the strong beams 
propagate parallel and anti-parallel to the probe field. In this experimental set-up, we simply 
exchange the co-ordinates y and z in the expressions for the strong fields in Eq. ( 1II.5j) , gaining 
a corresponding y-axis strong- field Rayleigh length, y r and Gouy phase g)O = tan _1 (|//|/ ri o). 
To the probe field, we add defocusing terms inside the integral, which would allow us to 
consider the case w p>0 < wo,o- The diffracted field, Ed(rd) is then given by the following 
expression: 

I aE p 



E d (r d ) 



2VX + (V^ + V^), 



(111.20) 
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(111.21) 
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where we have also introduced a phase difference term A-^o = ^0,2 — i/>o,i between the 
absolute phases of the two strong beams, which turns out to have negligible effect for the 
same reasons as separating the beams in the longitudinal direction and is correspondingly 
set to zero. The vectors u[ and u' 34 are: 



( (1 + ^)cos# + ^(^cosfl + l^sintf) \ 



Ui 



r d / ■ r d \ r d ■ t r d 

4 r d r d r d 4 r d r d 

sm6 + ^ cos9 + |((^) 2 + 1) sinfl 



(111.23) 
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\ - I)sin0 + 2^cos0- fp) 2 sin# / 

\ 4 V r d / r d r d 4 V r d / / 



As in the previous case, we set (xd/?"d), (zd/?"d) - ► 0, (yd/r^) — )■ 1 in Eq. ( 1III.23I) and 
Eq. (1III.24[) which removes the latter vector completely, eliminating any contribution from 
Eo.2( r ? £)) the strong-field beam with wavevector parallel to that of the probe. This result is 
expected from the general property of a plane wave that it does not polarise the vacuum. In 
this geometry, we obtain for the polarisation ifj and the ellipticity e the following expressions: 
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n sin 2^ / f Vf Vf 



ISA 2 



3/r,p 
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(111.25) 
(111.26) 

We can compare these to existing results arrived at by Heinzl et al. [39| when we take 
a = b = and the two limits: the refractive-index and the crossed-field limit. The first is 
obtained when we take y& — > (near region), in a regime where if) becomes linear with y& 
and therefore disappears, and e converges to a constant. The crossed-field limit corresponds 
to a constant strong field, i.e. uon — > 0, which we can achieve when we let the counter- 
propagating pulse be e.g. of the form of a cosine. This ensures that neither the strong 
electric nor magnetic field disappears in this limit, so that we can keep the normalisation 
used in Eqs. flIII.25j) and (IIII.26I) . To be consistent, the time-averaging procedure which 
removes evanescent waves must be repeated with the precondition that ojq = 0. Then Eq. 
(1III.26D tends to the result in [391]: 
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The only difference to the formula in 39] is that we have incorporated the focusing of the 
strong- and probe- fields, which automatically generates the effective interaction length yo 
of the beams. 



Another feature which is different here, is that we allow the strong-field wave to be positioned 
off-axis. We showed and explained how this increases the ellipticity and polarisation in the 
double-slit set-up, and in this single-shaft experiment with just one beam, one acquires a 
similar result (see Fig. [8]). For the same experimental parameters as in Fig. (T5]) but with 
yd = 50 cm and a/w 0j o = 10, we achieve a modest increase in the ellipticity of 1.3 over single 
strong-beam values. We mention here, that one could also form a double-shaft geometry 
which leads to the same relative increase as for the off-axis single-shaft one. As Xd = Zd = 0, 
this can be understood as a result of the symmetry of the set-up, which can also be seen in 
Eq. ( lIII.21j) . being symmetric in x and a in this limit. 



IV. CONCLUSION 



A main focus of this paper was to extend the results derived in [40( to incorporate more 
features applicable to experiment. One development has been to extend into the far- field 
region, the range in which polarisation rotation and ellipticity formulae are valid. These 
results were calculated for two different geometries: double-slit and single/double-shaft. 
Another addition has been to include a separation of the strong-field beams. This non-trivial 
beam geometry in conjunction with higher-order Fresnel diffraction terms was shown to 
increase polarisation rotation and ellipticity values for a range of beam parameters in the 
double-slit case by a factor of 1.4, and in the single/double-shaft case by a factor of 1.3. 
Although these increases are relative to the values at zero beam separation, we acknowledge 
that since the overall accuracy of the calculation is ~ l/V, some care should be taken 
in interpreting these results. By calculating the diffraction pattern resulting from the 
interference between the polarised vacuum and probe signals, we have illuminated another 
possible route to measuring laser-induced VPEs. For experimental parameters comfortably 
attainable at the upcoming X-FEL and ELI facilities, we have shown how approximately 
10 -5 of the incident photons can be diffracted, with around two photons per shot of the 
lasers being diffracted into regions where the vacuum signal is higher than the probe 
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background. However, we stress that an increase would also be observed using ELI with 



a table-top X-ray laser such as e.g. in 



54| where a beam of frequency 29 nm was used. 



These, in principle measurable diffraction vacuum polarisation effects, would be the first 
evidence of non-linear vacuum polarisation in laser fields. 
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Appendix 



The volume integral from Eq. (1II.15j) can be integrated in the x and y co-ordinates to give: 
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and included all four integrals with: 



< 



1 if k = 1, 
-1 if jfe = 2, 



if fc = 3,4, 



and 



ft 



if fc = 1,2, 

1 if k = 3, 
-1 if As = 4. 
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FIG. 1: The double-slit experimental setup. A monochromatic Gaussian probe beam with electric 
and magnetic field vectors E p and B p respectively, wavevector k p , linearly polarised at an angle 
9 to the x-axis in the x-z plane and with a waist w P; o, much greater than the strong-field beam 
waist, wo,0) impinges and is perpendicular to two, parallel, counter-propagating, monochromatic 
and Gaussian strong- field waves with amplitudes Eq / y/2 3> E p , electric and magnetic fields in the 
x-y plane Eo,i,Eo,2 and Bo,i,Bo,2 respectively, wavevectors ±ko = (0, 0, T^o), w o oj p , with foci 
at (x,y) = (a,b) and (x, y) = (—a, —b). The results of this process are then measured a distance 
r<i from the centre of the interaction region. 
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FIG. 2: The one-dimensional diffracted field along the x-axis that is predicted by first-order QED 
theory for the parameters: a/wo,o = 6,fo/wo,o = 0,wo.o = Ao = 0.8 ^m, A p = 0.4nm,w Pj o = 
100 /mi, 9 = it/2, y d = 1 m, 1$ = 10 24 Wcm -2 , is plotted with a solid line. The dashed line indicates 
the result obtained by using the simplified analytical approach based on Eq. ()III.3|) . 
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FIG. 3: The one-dimensional diffracted field along the z-axis i.e. along the axis of propagation 
of the strong-field beams, predicted by first-order QED theory is plotted with a solid line for the 
same physical parameters as in Fig. [2] but now with a/woo = 0. The dotted line represents the 
simplified analytical approach keeping quadratic terms in the exponential, described in the text. 
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FIG. 4: Plot a is of the quantity I p d, the vacuum-probe cross-term in time-averaged total intensity, 
at a distance « y& = 1 m, with the same experimental parameters as in the example of diffraction 
along the x-axis apart from a/wo,o = 12. Plot b is the logarithm to base 10 of this divided by the 
time-averaged probe intensity i.e. the logarithm of the signal-to-noise ratio for the same parameters. 
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FIG. 5: For a probe beam focused to w Pi o = 3.6 fim, w p (j/d) = 5.0 pm and a large enough detector 
hole radius p, the photon count from the vacuum-probe cross term iV p d(p) (solid-line) becomes 
comparable to that from the probe N p (p) = iV p (0) exp(— 2(p/w p (yd)) 2 ) (dashed-line) (here around 
2 diffracted photons) and greater than statistical noise from the probe ^N p (p) (dotted-line). 
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FIG. 6: The absolute value of the new polarisation and ellipticity expressions derived with probe 
defocusing terms and separated strong field beams as a function of the observation distance y^, 
for the same parameters as in Fig. ([2]) but with a/wo,o = 12. The dashed lines are the former 
analytical polarisation and ellipticity formulae for a = without probe defocusing terms, given in 



40 ], with the solid lines the new polarisation and ellipticity presented in this paper. Darker lines 



are drawn for polarisation rotation and lighter ones for ellipticity 
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FIG. 7: The ratios of \ip(a)/ip(0)\ (continuous line) and \e(a)/e(0)\ (dashed line) in the double-slit 
set-up for the same parameters as in Fig. ([2]) but with w Pj q = 200 fim. The ratio \ip(a)/ip(0)\ was 



chosen in preference to the ratio of |^( a )| with that in the original paper 
focusing terms doesn't make the latter a viable comparison. 



401 ] as introduction of 
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FIG. 8: The ratios of \ip(a)/ip(0)\ (continuous line) and |e(a)/e(0)| (dashed line) in the single-shaft 
set-up for the same experimental parameters as in Fig. ([2]) but with = 50 cm. 
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